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Abstract 

Isolated populations of Peromyscus truei comanche occur in the Tule and Palo Duro Canyon 
region of the Texas Panhandle. Four plausible routes that may explain the current distribution 
of the pinon mouse were investigated herein by comparing the phylogenetic relationships of P. t. 
comanche with eight other subspecies of P. truei from throughout the western and southwestern 
United States and the Baja Peninsula of Mexico. To determine the origin and affiliation of P. t. 
comanche populations, DNA sequence data were obtained from the mitochondrial cytochrome-# 
gene (1,143 bp) and displacement loop (971 bp) and analyzed in a phylogenetic context. Results 
of the Bayesian inference analyses indicated that samples of P. t. comanche genetically were undif¬ 
ferentiated from two of the eight subspecies of P. truei (P. t. nevadensis and P. t. truei). Overall, 
genetic divergence values, obtained from the cytochrome-# gene, ranged from 1.16% to 4.99% 
between P. t. comanche and the other eight subspecies of P. truei. Peromyscus truei comanche 
differed genetically from the nearest population ofP. t. truei (Glenrio, Texas) by 1.73%; however, 
samples off. t. comanche were differentiated from individual populations of P. t. truei from Mills 
Canyon, New Mexico, Black Mesa, Oklahoma, and Guadalupe Mountains National Park, Texas 
by 0.79%, 1.00%, and 2.00%, respectively. Although results of the phylogenetic analyses were 
unable to definitively explain the origin of P. t. comanche , levels of genetic divergence suggest 
a recent evolutionary history with samples from Mills Canyon, New Mexico. Results from the 
divergence dating analyses suggest that P. t. comanche from the Texas Panhandle diverged from 
a source-stock population in northeastern New Mexico approximately 0.71 mya. This infers 
an expansion route (Route II) from Mills Canyon, New Mexico, along canyonlands associated 
with the Canadian River system to the northeastern edge of the Llano Estacado, and subsequent 
isolation and divergence of P. t. comanche from its P. t. truei- like ancestor; whereas the closest 
extant populations of P. t. truei in Glenrio, Texas genetically were more divergent (1.73%), and 
also would have required dispersal across the inhospitable landscape (for P. truei ) of the Llano 
Estacado. Finally, P. t. comanche did not form a monophyletic group (to the exclusion of other 
samples of P. t. truei ) in the genetic analyses, however the significant levels of morphological 
differentiation reported in other studies justifies its retention as a valid subspecies. 
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Introduction 


Peromyscus truei occurs throughout the south¬ 
western United States and portions of Sonora and Baja 
California del Norte, Mexico (Hotfmeister 1981; Durish 
et al. 2004). Described by Blair (1943), P. t. coman¬ 
che geographically is isolated along the rocky, juniper 
slopes of northern and eastern peripheral edges of the 
Llano Estacado, specifically in Randall, Armstrong, 
and Briscoe counties of the Tule and Palo Duro Canyon 
areas of the Texas Panhandle (Blair 1943; Choate 1991, 
1997). P. t. comanche prefers more precipitous rocky 
areas along the escarpment and thus a microhabitat 
separation exists between P. t. comanche and other 
sympatric species such as P. attwateri and P. leucopus 
(Choate 1991, 1997; Choate et al. 1991). This region, 
referred to as the ‘Break of the Plains’, contains cedar 
forests on the canyon slopes and serves as a suitable 
habitat for a specialist rodent like the pinon mouse 
(Schmidly 1973). 

The occurrence and explanation for the distri¬ 
bution of P. t. comanche has proven to be enigmatic. 
Although several populations of P. t. comanche occur 
along the northeastern side of the Llano Estacado, 
(Choate 1991, 1997; Durish et al. 2004; Schmidly and 
Bradley 2016), the nearest additional populations of P. 
truei , excluding the Tule and Palo Duro Canyon locali¬ 
ties, are samples located in Deaf Smith County (south 
of Glenrio, Texas; Choate 1997; Choate et al. 1991) and 
Quay County (SE of Tucumcari, New Mexico; Findley 
et al. 1975; Hall 1981) at the interface of the western 
edge of the Llano Estacado and Canadian River Canyon 
and extreme western Oklahoma (Tesquite Canyon and 
Kenton region; Caire et al. 1989; Hall 1981). 

The aforementioned samples from eastern New 
Mexico, western Oklahoma, and northwestern Texas 
have been assigned to P. t. truei (Hall 1981; Choate 
1997; Choate etal. 1991; Schmidly and Bradley 2016). 
These localities are 120-180 km W and 300 kmNNW, 
respectively from the type locality of P. t. comanche 
near Tule and Palo Duro Canyons, Briscoe County, 
Texas. How P. t. comanche populated the northern 
and eastern edge of the Llano Estacado is unknown; 


however, given its proximity to populations of P. t. truei 
near the New Mexico/Texas border (circa Glenrio and 
Tucumcari) and western Oklahoma (circa Black Mesa), 
it is possible either the western or northwestern popu¬ 
lations of P. t. truei may have served as the ancestral 
population of P. t. comanche. 

Blair (1943) described P. t. comanche as a distinct 
species (P. comanche) based on a series of 92 specimens 
collected from the Tule and Palo Duro Canyon region 
of the Texas Panhandle. Blair’s description was based 
on the observation thatP comanche possessed: smaller 
auditory bullae; a shorter ear relative to hindfoot; and 
a significantly longer tail than did other members of 
the pinon deermouse, P. truei (Blair 1943). In Hoff- 
meister’s (1951) review of the P. truei species group, P. 
comanche was relegated to a subspecies of P. nasutus 
based on coloration, a uniquely shaped interparietal, 
and larger auditory bullae. Hotfmeister and de la Torre 
(1961) split P. difficilis andP nasutus into different spe¬ 
cies with P. n. comanche relegated to subspecific status 
within P. difficilis due to its relatively small ears, shorter 
tail, and small auditory bullae. Johnson and Packard 
(1974), based on results of an allozyme study, deter¬ 
mined that P. comanche genetically was most closely 
affiliated with P. truei. Further, although they found 
no apparent fixed allozymic differences between P. 
comanche and P. truei , they concluded that P. comanche 
should be considered as a species distinct fromP. truei. 
Conversely, based on morphologic and karyotypic data, 
Schmidly (1973) and Modi and Lee (1984) retained P. 
t. comanche as a subspecies of P. truei. 

Using haplotype patterns developed from restric¬ 
tion digests of the mitochondrial genome, DeWalt et 
al. (1993) argued that P. t. comanche should remain a 
subspecies of P. truei. Similarly, Durish et al. (2004) 
sequenced the mitochondrial cytochrome-# ( Cytb ) gene 
and reported that individuals ofP t. comanche differed 
from individuals of P. t. truei collected from New 
Mexico, Arizona, and California by a genetic diver¬ 
gence of 1.20%. Although Durish et al. (2004) agreed 
with Modi and Lee (1984), Dewalt et al. (1993), and 
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Schmidly (1973) that P. t. comanche should continue 
to be recognized as a member of P. truei, unfortunately, 
they were unable to examine specimens of P. t. truei 
from nearby populations in western Texas, eastern New 
Mexico, or extreme western Oklahoma to help establish 
a phylogenetic connection of P. t. comanche to the more 
western and northern populations. Furthermore, they 
did not offer any comments on whether P. t. comanche 
was a valid subspecies of P. truei or whether it should 
be subsumed into P. t. truei , the geographically closest 
population. 

Recent evidence presented by Rogers et al. (2019) 
suggested that the systematics of P. truei are more 
complex than reported in these previous studies. Their 
work, using DNA sequence and morphologic data, 
indicated that P. truei, as currently recognized, may 
be comprised of multiple species, including western 
and eastern phylogroups. Furthermore, based on their 
data, P. t. comanche was aligned with the more eastern 
phylogroup, which includes portions of P. t. truei and 
P. t. nevadensis (Rogers et al. 2019). 

In addition to the uncertain phylogenetic place¬ 
ment of P. t. comanche, the origin of the Tule and Palo 
Duro Canyon population is not clear, although there 
seem to be two possibilities, one of which is that P. 
t. comanche populations successfully dispersed via a 
“stepping-stone model”, presumably from nearby popu¬ 
lations in eastern New Mexico/western Texas or from 
western Oklahoma. Under this model, adequate habitat 
patches would facilitate movement from one population 
to another; thereby, successfully providing temporary 
refugia for subsequent dispersal and expansion to more 
optimal habitat (Yang et al. 2016). Kimura (1953) and 
Kimura and Weiss (1964) suggested that the stepping- 
stone model would explain movement of individuals 
along fragmented habitat and this scenario would sup¬ 
port documented distributions for P. t. comanche along 
the escarpments and edges of the Llano Estacado and 
into the Tule and Palo Duro Canyon area. 

Two possible scenarios can be proposed to ex¬ 
plain the present-day distribution of P. t. comanche. 
Under the first scenario, individuals presumably may 
have dispersed (via stepping-stone model) from western 
Texas, northeastern New Mexico or western Oklahoma. 
These putative routes of dispersal (IA, IB, II, III, and 
IV), or some combination thereof, are shown in Figure 


1. Route I (Glenrio, Texas to the Tule and Palo Duro 
Canyon region) represents the most plausible disper¬ 
sal route, although it would require individuals of P. 
t. truei to traverse 120 km across unsuitable habitat. 
Routes IIA and IIB (northeastern New Mexico to the 
Tule and Palo Duro Canyon region) depict potential 
pathways, with P. t. truei presumably dispersing either: 
IIA) northeast following the edges of the Canadian 
River system and, subsequently, south following the 
eastern edges of the Llano Estacado into Palo Duro 
Canyon; or IIB) following the Canadian River system 
to other tributary systems (major: Palo Duro and Tierra 
Blanca creeks; minor: Timber, Sunday, and Canyon 
Cita creeks; see Gould 1907) in a southeasterly direc¬ 
tion, directly leading into Palo Duro Canyon. Route 
III (Black Mesa to the Tule and Palo Duro Canyon 
region) depicts an alternative pathway of dispersal 
along the edges of the Wichita and Amarillo Mountain 
Uplifts to the northeastern edges of the Llano Estacado 
southward into Tule and Palo Duro Canyon. Route IV 
(Guadalupe Mountains National Park to the Tule and 
Palo Duro Canyon region) describes a final pathway 
along the southwestern and southeastern edges of the 
Llano Estacado northward to the southern distribution 
ofP t. comanche populations. 

A second scenario (vicarience) would suggest that 
P. t. comanche was a relictual population that originally 
was a widespread population and subsequently became 
isolated during the development of the Llano Estacado. 
Uplifts of the Rocky Mountains during the Cenozoic 
(70 million years ago [mya]), produced the materials 
underlying the Llano Estacado, and subsequent ero¬ 
sion and weathering to the Rocky Mountains by the 
Pecos and Canadian rivers generated the canyonland 
landscape in which P. t. comanche resides (Choate 
1991,1997). Initial separation due to these aforemen¬ 
tioned geographical barriers to dispersal along with P. 
truei specialist tendencies (e.g. preference for rocky 
outcroppings and canyon slopes), in recent times, may 
have further isolated populations of P. truei along the 
northeastern edge of the Llano Estacado followed by 
a divergence into P. t. comanche. 

To discern among these two scenarios, a combi¬ 
nation of mitochondrial markers (Cytb; displacement 
loop, D-loop) and phylogenetic analyses were used to 
determine the origin of P. t. comanche relative to other 
geographically proximal populations. The objectives 


4 


Occasional Papers, Museum of Texas Tech University 



Figure 1. Sampling localities for selected individuals of P. truei examined in this study. 
Samples from GenBank, including Rogers et al. (2019), are not depicted. Numbers contained 
in white circles indicate localities from Texas and surrounding states. Roman numerals (I, 
IIA, IIB, III, and IV) indicate hypothetical pathways for P. truei dispersal onto the Llano 
Estacado and surrounding canyonlands. Dashed lines indicate the distribution of P. t. truei 
and P. t. comanche. Refer to the Appendix for a complete list of specimens examined. 


of this study were to: 1) investigate the potential geo¬ 
graphic impact of the Llano Estacado on the current 
distribution of P. truei ,; and 2) determine if the Tule 
and Palo Duro Canyon populations are more closely 
aligned with populations to the west, northwest, or 


southwest, or are they a relict of a once more broadly 
distributed P. truei. In addition, molecular data were 
used to determine which of the scenarios, stepping- 
stone or vicariance, explains the disjunct distribution 
ofP. t. comanche. 
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Materials and Methods 


Sampling. —Individuals ( n = 162) represent¬ 
ing nine subspecies of P. truei were collected from 
naturally-occurring populations in the western United 
States and Mexico, 28 of which were borrowed from 
the Natural Science Research Laboratory, Museum of 
Texas Tech University and 134 were incorporated from 
other molecular studies (e.g., GenBank). Specimens 
were collected following methods outlined in the guide¬ 
lines of the American Society of Mammalogists (Sikes 
et al. 2016) and approved by the Texas Tech University 
Animal Care and Use Committee. To determine the 
source-stock of P. t. comanche , samples of P. t. truei 
from Mills Canyon, New Mexico; Black Mesa State 
Park and Nature Preserve, Oklahoma; Glenrio, Texas; 
and Guadalupe Mountains National Park, Texas, were 
used to represent the closest known populations to 
the current distribution of P. t. comanche. For the 
Cytb dataset (see Appendix), reference sequences ( P .! 
attwateri , gratus, nasutus, and truei) were obtained 
from the following studies: 117 individuals (Rogers 
et al. 2019), 9 individuals (Durish et al. 2004), 3 indi¬ 
viduals (Rodhouse et al. 2010), 3 individuals includ¬ 
ing P. attwateri as outgroup taxon (Tiemann-Boege 
et al. 2000), 1 individual (Bradley et al. 2004), and 1 
individual (Smith and Patton 1999). For the D-loop 
dataset, 35 individuals (see Appendix) were sequenced 
following the methods of Castro-Campillo et al. (1999). 
All sequences ( Cytb and D-loop) generated herein were 
deposited in GenBank and are listed in the Appendix. 

DNA sequencing. —For Cytb sequences, genomic 
DNA was isolated from O.lg of frozen liver tissue 
using the Qiagen DNeasy kit (Qiagen Inc., Valencia, 
California). For all samples the full length Cytb gene 
(1,143 bp) was amplified using polymerase chain reac¬ 
tion (PCR) methodology following Saiki et al. (1988) 
with the amplification primers: LGL765 (forward, 
Bickham et al. 1995) and LGL766 (reverse, Bickham 
et al. 2004). HotStarTaq (Qiagen Inc., Valencia, Cali¬ 
fornia) was used in the PCR along with the following 
thermal profile: hot start of 80°C, initial denaturation at 
95°C for 2 min, followed by 34 cycles of denaturation 
at 95°C for 30 s, annealing (range: 44-45°C) for 45 s, 
extension at 73°C for 1 min, and a final extension at 
73°C for 15 min. 

Tissue samples were unavailable for the speci¬ 
mens from Glenrio and Guadalupe Mountains National 


Park, Texas; therefore, a small piece of a toe (approxi¬ 
mately 5 mg excluding the toenail) was obtained from 
four (two from each locality) museum voucher speci¬ 
mens. DNA was isolated following a modification of 
methods reported in Curry and Derr (2019). Toes 
initially were cleaned using a 95% ethanol rinse and 
then immediately treated with UV irradiation for 5 
min. Each toe clip was washed with ddH^O and incu¬ 
bated at 56°C for 15 minutes (3 repetitions). Genomic 
DNA was isolated using the Quick-DNA Universal Kit 
(Zymo Research, Irvine, California). A 423 bp frag¬ 
ment (position 400-823 bp aligned) of the Cytb gene 
was amplified using PCR methods with primers 400F 
(Edwards et al. 2001) and 700H (Peppers and Bradley 
2000), HotStarTaq (Qiagen Inc., Valencia, California) 
and the same thermal profile as described above for the 
tissue samples with the exception that the annealing 
temperature was 50°C. 

All PCR products were purified with either 
ExoSAP-IT or a Qiagen PCR Purification Kit (Ap¬ 
plied Biosystems, Foster City, California and Qiagen 
Inc., Valencia, California). Cycle sequence primers 
included LGL765, LGL766, 870R, and FI (Peppers 
et al. 2002; Whiting et al. 2003) and subsequent cycle 
sequencing reactions were purified using the ABI Prism 
Big Dye version 3 terminator ready reaction mix (Ap¬ 
plied Biosystems, Foster City, California). Sequenc¬ 
ing reactions were purified using Sephadex columns 
(Princeton Separation, Adelphia, New Jersey) and cen¬ 
trifugation, followed by dehydration and resuspension 
in formamide. Purified products were analyzed on an 
ABI Prism 310 automated sequencer (Biotechnology 
Resource Center, Institute of Biotechnology, Cornell 
University, Ithaca, New York). Resulting sequences 
were aligned and proofed with Sequencher 4.10.1 soft¬ 
ware (Gene Codes Corporation, Ann Arbor, Michigan) 
and chromatograms were inspected to authenticate any 
base changes. 

The mtDNA D-loop was amplified using a similar 
methodology as described above. The differences in 
sequencing methods for D-loop are described below. 
Primers utilized to amplify the full-length D-loop (971 
bp) were 2340-4 (forward, Bickham et al. 1995) and 
2340-5 (reverse, Castro-Campillo et al. 1999). Ther¬ 
mal profiles for PCR were as follows: a hot start of 
80°C, initial denaturation at 95°C for 2 min, followed 
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by 35 cycles of denaturation at 95°C for 30 s, anneal¬ 
ing (range: 48-49°C) for 45 s, and extension at 72°C 
for 1 min, with a final extension at 72°C for 15 min. 
Primers used to cycle sequence the products included 
2340-4, 2340-5, 500F, and 1115 (Mendez-Harclerode 
et al. 2005). 

Phylogenetic analyses. —For the Cytb gene, a par¬ 
simony analysis was conducted using PAUP* 4.0al66 
(Swofford 2003) on an initial dataset containing 523 
individuals. Nucleotide positions were treated as un¬ 
ordered, discrete characters with four possible states: 
A, C, G, and T. Verification of taxonomic assignment, 
elimination of duplicate haplotypes, and confirmation 
of monophyletic clades resulted in a final dataset of 
162 individuals for subsequent analyses. Peromyscus 
attwateri, the sister taxon to P. truei (Durish et al. 2004), 
was incorporated as the outgroup taxon. 

jModelTest-2.1.10 (Darriba et al. 2012) identified 
GTR+G and TIM3+I+G as the most appropriate model 
of evolution for Cytb and D-loop, respectively. To 
perform likelihood analysis under a Bayesian inference 
model, MrBayes v.3.2.6 (Ronquist et al. 2012) was run 
with the following parameters: 2 independent runs with 
four chains, one cold and three heated (MCMCMC), 
10 million generations, and sample frequency every 
1,000 th generation from the last nine million generated. 
A consensus tree (50% majority rule) was constructed 
from the remaining trees, and posterior probability 
values (PPV) were estimated to provide nodal support. 
PPV >0.95 were considered indicative of nodal sup¬ 
port. Genetic distance values were calculated using the 
Kimura 2-parameter model of evolution (Kimura 1980) 
to assess levels of divergence between populations ofP 
truei (Arizona, Baja California, California, Colorado, 
Idaho, New Mexico, Oklahoma, Oregon, Texas, and 
Utah) and P. t. comanche (Texas) following Baker and 
Bradley (2006). 

Genetic divergence and divergence dating. —Two 
methods were used to examine genetic divergence and 


Phylogenetic analyses. —Complete nucleotide 
sequences from the mitochondrial Cytb gene (1,143 bp) 
and D-loop (971 bp) were obtained from 162 and 35 


molecular dating. In the first method, the Kimura 2-pa¬ 
rameter model of evolution (Kimura 1980) was utilized 
to estimate genetic distances among selected taxa and 
haplotypes. The resulting values were used to examine 
levels of genetic divergence pertaining to the genetic 
species concept outlined in Bradley and Baker (2001) 
and Baker and Bradley (2006). A second method (mo¬ 
lecular clock test; MEGA v. 10, Kumar et al. 2018) was 
employed to define the appropriate molecular timeline 
and to determine whether to accept or reject a stringent 
molecular clock. Divergence dates for subspecies ofP. 
truei were estimated from the Cytb dataset (obtained 
in this study and GenBank, see Appendix) using the 
program BEAST v2.6.1 (Bouckaert et al. 2014). Di¬ 
vergence date estimates were not determined for the 
D-loop dataset due to lack of in-depth geographic and 
taxonomic sampling. 

Fossil calibrations were placed on the P. truei 
node, based on the fossil date (~2.7 mya) from the most 
recent common ancestor to P. attwateri (Dalquest 1962; 
Karow et al. 1999) following methods outlined in previ¬ 
ous studies (Schenk et al. 2013; Bouckaert et al. 2014; 
Ordonez-Garza et al. 2014; Platt et al. 2015; Sullivan 
et al. 2017; Bradley et al. 2019). A Yule tree prior was 
used for the BEAST analysis and a prior lognormal 
distribution was placed on root height to constrain the 
divergence date estimates of the overall tree to ~2.7 
mya with a o value of 0.5 and to reflect the uncertainty 
in the fossil record. Outgroup taxa used were the same 
as the Bayesian analyses. Optimization of the analy¬ 
sis and final parameters were evaluated using initial 
test runs with the following parameters: GTR+I+G 
(nucleotide substitution model), 1.0 x 10 7 generations, 
and 10% burn-in. Two final runs of 2.0 X 10 7 genera¬ 
tions were analyzed with log and tree files, which were 
then combined to generate divergence date estimates 
and produce a maximum clade credibility tree. Tracer 
and TreeAnnotator (Bouckaert et al. 2014) was used 
to examine for sufficient mixing, convergence stability, 
and effective sample size >200 for all parameters and 
to estimate the optimal phylogenetic tree, respectively. 


individuals, respectively. An initial Bayesian analysis 
used P. attwateri as the outgroup taxon in addition to 
P. difficilis , P. nasutus, and P. gratus. A secondary 


Wright et al.—On the Origin of Peromyscus truei comanche 


1 


Bayesian analysis only used P. gratus (sister taxon) as 
the outgroup. Both analyses produced similar topolo¬ 
gies but only the first is reported herein. The topology 
of the phylogenetic tree obtained from the Bayesian 
analysis of the Cytb dataset indicated nodal support 
(PPV > 0.95) primarily for basal nodes; consequently, 
unsupported clades were collapsed so that only major 
supported clades were depicted (Fig. 2). Members of 
each species (P. difficilis , P nasutus, P. gratus , and P. 
truei ) were arranged as independent, monophyletic 
clades. Within P. truei , two supported clades were 
identified that corresponded to the western haplotype 
(Arizona, Baja California, California, Idaho, Nevada, 
Oregon, and Utah) and eastern haplotype (Arizona, 
Colorado, Oklahoma, New Mexico, Texas, and Utah) 
similar to that depicted in Rogers et al. (2019). How¬ 
ever, relationships between members of P. truei sub¬ 
species and P. t. comanche revealed little to no genetic 
differentiation and no evident population structure. The 
Bayesian phylogenetic tree, generated using D-loop 
DNA sequences (Fig. 3), produced a topology in which 
18 of the 30 nodes were supported with Bayesian PPV 
> 0.95. Similar to the Cytb phylogeny, members of 
each species group were arranged into four supported 
clades and within P. truei there was little to no genetic 
differentiation and structure. 

Genetic distances. —For selected taxa, average 
genetic distance values (Table 1) were obtained from 
the Kimura 2-parameter model of evolution (Kimura 
1980). Interspecific comparisons ranged from 7.21% 
to 13.30% {Cytb) and 5.58% to 11.10% (D-loop) 
whereas intraspecific comparisons for both Cytb and 
D-loop ranged from 1.00% to 3.00%. Average genetic 
distances of Cytb between P. t. comanche and the other 
eight subspecies of P. truei ranged from 1.16% (P t. 


Results of the Cytb analyses indicated two sup¬ 
ported clades representing a geographical designation 
(eastern and western) conforming to that discussed by 
Rogers et al. (2019). All individuals ofP. t. comanche 
were included in the eastern clade along with some rep¬ 
resentatives of P. t. nevadensis and P. t. truei. However, 
other samples of P. t. nevadensis and P. t. truei were 
included as members of the western clade along with P. 
t. chlorus, P. t. gilberti, P. t. lagunae, P. t. martirensis, P. 


truei east) to 4.99% {P. t. chlorus). Subspecific com¬ 
parisons for the four focal populations of P. t. truei 
(Mills Canyon, Black Mesa, Glenrio, and Guadalupe 
Mountains National Park) to populations ofP t. coman¬ 
che (Tule and Palo Duro Canyon region) ranged from 

O. 79% to 2.00% (Cytb); subspecific comparisons for 
D-loop were only available between P. t. comanche and 

P. t. truei (Mills Canyon = 1.71%) and P. t. comanche 
and P. t. truei (Black Mesa = 1.81%). 

Divergence dating. —Molecular clock tests (Ku¬ 
mar et al. 2018) rejected a strict molecular clock for the 
Cytb dataset; therefore, a relaxed log normal molecular 
clock was used in subsequent analyses. Mean rates of 
evolution (p < 0.05) between clades depicted in the 
Bayesian analysis yielded a mean rate of 0.3 substitu¬ 
tions per site per million years (95% highest posterior 
density [HPD]: 0.66 to 5.20) for Cytb (Fig. 4). The Yule 
birth rate was estimated to be 4.25 (95% HPD: 0.98 to 
8.83). The divergence dating analysis indicated that the 
initial radiation of the western and eastern clades of P. 
truei was approximately 1.29 mya. The divergence date 
estimates obtained from the Cytb dataset indicated that 
a P. truei-\ike ancestor diverged from a P. gratus-Vke, 
ancestor approximately 1.85 mya. Within P. truei, the 
western and eastern clades appear to have diverged 1.29 
mya, with the western haplotype corresponding to eight 
subspecies (P. t. chlorus, P. t. gilberti, P. t. lagunae, P. 
t. martirensis, P. t. montipinoris, P. t. nevadensis, P. t. 
preblei, and P. t. truei) and the eastern haplotype cor¬ 
responding to three subspecies {P. t. comanche, P. t. 
nevadensis , and P. t. truei). Relative to this study, the 
divergence of the three closely related subspecies {P. 
t. comanche, P. t. nevadensis , and P. t. truei) contained 
in the eastern clade occurred approximately 0.81 mya. 


t. montipinoris , and P. t. preblei. Based on the results 
of Rogers et al. (2019) and data presented herein, in¬ 
dividuals composing the eastern clade probably should 
be recognized as a single taxonomic unit distinct from 
those in the western clade. DNA sequences from D- 
loop did not contribute substantial information to the 
phylogenetic relationships among the nine subspecies 
of P. truei due to either the low levels of genetic di¬ 
vergence or lack of taxonomic sampling; consequently, 
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-| P. gratus (n = 2) 


— P. difficilis (n = 1) 
-\ P. nasutus (n = 5) 


P. attwateri (n = 1) 


P. t. truei (n = 43) 

P. t. nevadensis (n = 15) 
P. t. montipinoris (n = 16) 
P. t. gilberti (n = 16) 
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P. t. preblei (n = 1) 
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P. t. comanche (n = 9) 

P. t. nevadensis (n = 5) _ 
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Figure 2. Phylogenetic tree generated from cytochrome-/) sequences and Bayesian inference analyses (MrBayes 
3.2; Ronquist 2012) and the GTR+G model of evolution. Posterior probability values > 0.95 are indicated by 
an asterisk and depict nodal support. Designation of western and eastern haplotype are consistent with clades 
reported in Rogers et al. (2019). 
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Figure 3. Phylogenetic tree generated using displacement-loop sequences and Bayesian inference analyses 
(MrBayes 3.2; Ronquist 2012) and the TIM3+I+G model of evolution. Posterior probability values > 0.95 are 
indicated by an asterisk and depict nodal support. 
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Table 1. Average genetic distances estimated using the Kimura 2-parameter model of evolution 
(Kimura 1980) reported for select taxonomic groups in this study. East and west designations refer 
to the haplotypes defined by Rogers et al. (2019). Populations assigned to various subspecies and 
localities are depicted in the Appendix. Guadalupe Mountains National Park was abbreviated to 
GUMO. 


Peromyscus comparsion 

Cytb 

D-loop 

Between species 

P. attwateri versus P. difficilis 

8.50 

7.37 

P attwateri versus P. nasutus 

8.01 

8.02 

P attwateri versus P. gratus 

11.40 

10.50 

P. attwateri versus P. truei 

13.30 

9.49 

P. attwateri versus ingroup 

13.10 

9.27 

P. difficilis versus P nasutus 

7.21 

5.58 

P. difficilis versus P gratus 

10.80 

9.86 

P difficilis versus P. truei 

13.20 

10.90 

P. difficilis versus ingroup 

13.00 

10.10 

P. nasutus versus P. gratus 

11.50 

10.60 

P. nasutus versus P truei 

12.00 

11.10 

P. nasutus versus ingroup 

12.00 

11.10 

P. gratus versus P. truei 

10.40 

8.19 

Within species 

P. attwateri 

- 

- 

P. difficilis 

- 

- 

P. nasutus 

1.00 

1.00 

P gratus 

2.00 

3.00 

P. truei 

3.00 

2.00 

P. truei (west) 

1.00 

- 

P. truei (east) 

1.00 

- 

Between subspecies 

P t. comanche versus P. t. montipinoris 

4.81 

- 

P. t. comanche versus P t. gilberti 

4.89 

- 

P. t. comanche versus P. t. chlorus 

4.99 

- 

P. t. comanche versus P t. martirensis 

4.58 

- 

P. t. comanche versus P t. lagunae 

4.79 

- 

P t. comanche versus P. t. preblei 

4.88 

- 

P. t. comanche versus P t. nevadensis 

3.99 

- 

P. t. comanche versus P. t. nevadensis (west) 

4.88 

- 

P. t. comanche versus P t. nevadensis (east) 

1.32 

- 

P. t. comanche versus P. t. truei 

3.79 

- 
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Table 1. (cont.) 


Peromyscus comparsion 

Cytb 

D-loop 

P. t. comanche versus P. t. truei (west) 

4.82 

- 

P. t. comanche versus P. t. truei (east) 

1.16 

- 

P. t. truei (west) versus P t. truei (east) 

4.79 

- 

Within subspecies 

P. t. comanche 

1.00 

1.00 

P. t. truei 

3.00 

- 

P t. truei (west) 

1.00 

- 

P. t. truei (east) 

1.00 

- 

Between populations 

P. t. comanche versus P. t. truei (Glenrio, TX) 

1.73 

- 

P. t. comanche versus P t. truei (Black Mesa, OK) 

1.0 

1.81 

P. t. comanche versus P. t. truei (Mills Canyon, NM) 

0.79 

1.7 

P. t. comanche versus P t. truei (GUMO, TX) 

2.00 

- 



P. attwateri (n = 1) 


P. t. lagunae (n = 9) 

P. t. martirensis (n = 2) 

P. t. truei (n = 43) 

P. t. nevadensis (n = 15) 
P. t. montipinoris (n = 16) 
P. t. gilberti[n = 16) 

P. t. chlorus (n = 1) 

P. t. preblei[ n = 1) 


Western 

Haplotype 


P. t. truei (n = 36) 

P. t. comanche (n = 9) 
P. t. nevadensis (n = 5) 


Eastern 

Haplotype 


P gratus (n = 2) 


P. difficilis (n = 1) 


P. nasutus (n = 5) 


0.3 substitutions/site/mya 

To ' To 1 Jo 1 To 1 To ' Omya 
Figure 4. Time-calibrated phylogenetic tree modified from that depicted in Figure 2 with the superimposition of 
results from the BEAST analysis (Version 2.6.1, Bouckaert et al. 2014) using the mitochondrial cytochrome-/) 
gene dataset. Divergence date estimates are indicated along the x-axis in millions of years. Error bars (gray 
rectangles) represent the 95% highest posterior density for node height. 
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the majority of the discussion is based on conclusions 
from the Cytb dataset. 

The three populations (Fig. 2; Table 1) associ¬ 
ated with the Tule and Palo Duro Canyon region did 
not form a monophyletic assemblage; nor did they 
demonstrate geographical affinity to other populations 
of P. truei from northeastern New Mexico, western 
Oklahoma, or western Texas. However, genetic diver¬ 
gence values suggested a recent evolutionary history of 
haplotypes associated with P. t. comanche to P. t. truei 
haplotypes from Mills Canyon, New Mexico (0.79%) 
and Black Mesa, Oklahoma (1.00%). This result was 
unexpected given that the closest P. truei population to 
P. t. comanche was the P. t. truei samples from Glenrio, 
Texas (located 120 km to the west) which differed by 
1.73%. It is of interest that the Glenrio, Texas, popula¬ 
tion is by far the closest geographically extant popula¬ 
tion ofP. t. truei to those populations of P. t. comanche 
in the Tule and Palo Duro Canyon systems; however, 
a direct connection of Glenrio, Texas, to the Tule and 
Palo Duro Canyon systems would require dispersal 
across an intervening region of the High Plains that has 
been deemed inhospitable habitat for P. truei by Choate 
(1991, 1997). Further, results from divergence dating 
analyses suggested that haplotypes associated with 
samples of P. t. comanche from the Texas Panhandle 
possibly diverged from a source-stock population in 
northeastern New Mexico approximately 0.71 mya. 
This infers a possible expansion route (Fig. 5) from 
Mills Canyon, New Mexico, through canyonlands 
associated with the Canadian River system along the 
northeastern edge of the Llano Estacado, and subse¬ 
quent isolation and divergence of P. t. comanche from 
its P. t. truei-Uke ancestor. At this time, we cannot 
discern between routes IIA (along the canyonlands of 
the northeastern edge of the Llano Estacado) or IIB 
(along the southern major and minor tributaries of the 
Canadian River system). Although genetic divergence 
values associated with haplotype similarities between 
Black Mesa, Oklahoma, and Tule and Palo Duro 
Canyon were only slightly larger at 1.00% than those 
obtained from a comparison of haplotypes associated 
with P. t. truei (Mills Canyon, New Mexico, to the 
Tule and Palo Duro Canyon region) at 0.79%, route 
III (255 km) from Black Mesa, Oklahoma, to Tule and 
Palo Duro Canyon is unlikely due to absence of extant 
populations ofP truei along suitable habitat. Similarly, 
a relatively high level of genetic divergence (2.00%) 


between haplotypes associated with P. t. truei between 
Guadalupe Mountains National Park, Texas, and Tule 
and Palo Duro Canyon region in combination with a 
large intervening geographical distance between locali¬ 
ties (445 km) makes route IV unlikely. 

Divergence dating analyses indicated that all 
five focal localities diverged within the last 0.71 mya; 
therefore, it is reasonable to assume they were part of 
a continuous population across New Mexico, western 
Oklahoma, and western Texas. There is historical evi¬ 
dence for a mesic juniper-oak-pinyon pine woodland 
extending from the Rocky Mountains to the Texas/ 
Mexico border. However, in recent times (0.011 to 
0.008 mya), the lower elevation habitats transitioned 
to a more xeric juniper-dominated landscape, as shown 
by evidence from packrat middens and pollen (Hafsten 
1961; Wells 1966; Johnson and Packard 1974; Wells 
1977; Riskind and Van Devender 1979; Van Devender 
1990; Stangl etal. 1994; Abbott 1996; Bartlema2001). 
In examining the extant P. truei populations, Glenrio, 
Texas, Black Mesa, Oklahoma, and Palo Duro Canyon 
are more typical of the xeric juniper habitat whereas 
Mills Canyon, New Mexico, and Guadalupe Mountains 
National Park, Texas, populations are more similar 
to habitat associated with P. truei populations of the 
western distribution. Data generated herein do not 
show a modem-day connection (stepping-stone model) 
to Glenrio, Texas, and Black Mesa, Oklahoma, to the 
Tule and Palo Duro Canyon regions. Instead, the 
genetic divergence data suggests the most recent gene 
flow occurred from the Mills Canyon, New Mexico, 
population of P. t. truei to the Tule and Palo Duro Can¬ 
yon population representing P. t. comanche. Climatic 
changes producing a more xeric habitat then isolated the 
Tule and Palo Duro Canyon region population during 
the last 0.011 mya, resulting in a relictual population 
that eventually diverged morphologically to the extent 
that P. t. comanche was recognized as being distinct 
from P. t. truei. 

It may be necessary to incorporate advanced 
genomic techniques to resolve the phylogenetic 
relationships among the nine subspecies of P. truei and 
to describe the evolutionary history of P. t. comanche. 
Further, sampling from additional localities in eastern 
New Mexico may provide more information on the 
presence of the isolated meta-populations of P. t. 
comanche. Although there is no substantial genetic 
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Figure 5. Map depicting the four possible routes from putative source-stock populations of P. t. truei. 
The solid line indicates the most likely pathway (Route IIA or IIB) supported by the Kimura (1980) 
average genetic divergence values depicted in boxes. Dashed lines represent the three alternate routes 
(I, III, and IV) discussed in the text. Medium gray shading represents the presumed distribution of 
P. t. truei. Dark gray shading represents the distribution of the P. t. comanche population. Localities 
examined are abbreviated as follows: B = Black Mesa, Oklahoma; M = Mills Canyon, New Mexico; 
D = Glenrio, Texas; G = Guadalupe Mountains National Park, Texas; and P = Tule and Palo Duro 
Canyon, Texas. 
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differentiation between P. t. comanche and members of 
the eastern clade of P. t. truei (1.16%), and considering 
the isolated nature and morphological divergence as 
depicted by Blair (1943), Hoffmeister (1951), and 


Schmidly (1973), it seems prudent to continue to 
recognize P. t. comanche as a subspecies until additional 
data are available. 
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Appendix 

Specimens examined. —Specimens examined in this study are listed below. For each specimen, the general 
collecting locality is provided; all specimens were collected from the United States unless otherwise noted. For 
most taxa, the museum catalog number (abbreviations for museum acronyms follow Hafner et al. 1997) and 
GenBank accession numbers for Cytb , and D-loop are provided in parentheses and are separated by slashes, 
respectively. Multiple specimens from the same locality are separated by a semicolon. However, for samples 
of P. truei comanche and P. t. truei, the assigned locality, catalog number, and GenBank accession numbers for 
Cytb and D-loop are listed in parentheses and separated by slashes, respectively. Abbreviations are as follows: 
Carnegie Museum of Natural History (CM); Museum of Texas Tech University (TTU); Museum of Vertebrate 
Zoology at Berkeley (MVZ); Centro de Investigaciones Biologicas del Noroeste (CIB); University of Utah, Natural 
History Museum of Utah (UMNH); Brigham Young University, Monte L. Bean Life Science Museum (BYU); 
and University of Washington, Thomas Burke Memorial Washington State Museum (UWMB). The following 
specimens were sequenced for both Cytb and D-loop unless otherwise noted (NA). 

Peromyscus attwateri .—Oklahoma; McIntosh Co., 3.1 mi E Dustin (TTU55688/AF155384/N/A); Texas; Garza 
Co., 10 mi S Post (TTU36088/NA/AF081492). 

Peromyscus difficilis amplus. —MEXICO: Tlaxcala; 18 kmN, 9 kmE Apizaco (CM57108/AY387488/MN164782). 

Peromyscusgratusgentilis. —MEXICO: Durango; 3.8 mi WCoyotes, Hacienda Coyotes (TTU81621/AY322507/ 
MN164788). 

Peromyscus gratus. —New Mexico; Catron Co., Reserve, Gila National Forest (TTU138339/MN022892/ 
MN164789). 

Peromyscus nasutus. —New Mexico; Sandoval Co., 4 mi S, 3 mi W Bernalillo (TTU36097/MN022893/MN164784). 

Peromyscus nasutus. —Oklahoma; Cimarron Co., Black Mesa Nature Preserve, (TTU136503/MN022894/ 
MN164785; TTU138362/MN022895/MN164786). 

Peromyscus nasutus griseus. —New Mexico; Lincoln Co., 4 mi S Carrizozo (TTU78401/AF155399/MN164783). 

Peromyscus nasutus nasutus. —Texas; Jeff Davis Co., Mt. Livermore Preserve (TTU78316/AY376426/MN164787). 

Peromyscus truei chlorus. —California; San Bernardino Co., Cactus Flat, San Bernardino Mts. (MVZ 198708/ 
MK871875/NA). 

Peromyscus truei comanche. —Texas, Armstrong Co., 0.75 mi N, 6.25 mi E Wayside (Locality 3; TTU61541/ 
MN022908/MN164806; TTU61543/AY376428/MN164804; TTU61547/AY376429/MN164805); Briscoe Co., 
3 mi N Quitaque, Caprock Canyons (Locality 1; TTU47286/MN022915/MN164812; TTU47296/AY376430/ 
MN164808; TTU47300/MN022913/MN164810); 6 miN, 4 mi WSilverton (Locality 2; TTU36056/MN022912/ 
MN164809; TTU36060/MN022909/MN164815; TTU36064/MN022911/MN164814; TTU36068/MN022910/ 
MN164813); Caprock Canyons State Park (TTU74991/AY376431/MN164807; TTU69561/MN022914/ 
MN164811). 

Peromyscus trueigilberti .—California; Alameda Co., Strawberry Canyon (MVZ157329/AF108703/NA); Straw¬ 
berry Canyon, below Botanical Gardens (MVZ157330/MK871830/NA); Mariposa Co., 5.7 mi SE Coulterville 


18 


Occasional Papers, Museum of Texas Tech University 


(MVZ208171/MK871850/NA); Hunter Valley Mountain (MVZ208172/MK871851/NA; MVZ208173/MK871852/ 
NA; MVZ208176/MK871853/NA; MVZ208181/MK871857/NA); Blackstone Creek, 6.5 mi NE Coulterville 
(MVZ208184/MK871858/NA; MVZ208185/MK871859/NA; MVZ208187/MK871861/NA; MVZ208188/ 
MK871862/NA; MVZ208189/MK871863/NA); Monterey Co., Arroyo Seco, 7 mi SW Greenfield (MVZ195335/ 
MK871867/NA;MVZ195337;MK871868/NA;MVZ195338/MK871869/NA); ShirttailCanyon,4.8miE Soledad 
(MVZ195341/MK871870/NA). 


Peromyscus truei lagunae .—MEXICO: Baja California Sur; Valle de la Laguna, Sierra de la Laguna (CIB10956/ 
MK872276/NA; CIB 10959/MK872277/NA); Agua de San Antonio, 9 km N, 26 km E Todos Santos (CIB 10962/ 
MK872280/NA; CIB 10964/MK872281/NA; CIB 10965/MK872282); Los Pinitos, 17.5 km W Santiago (CIB 10966/ 
MK872283/NA); Palo Extrano, Sierra de la Laguna (CIB10969/MK872286/NA; CIB10972/MK872288/NA); 4 
km N, 22.5 km W Santiago (CIB10982/MK872293/NA). 


Peromyscus truei martirensis .—MEXICO: Baja California; Laguna Juarez (CIB3357/MK872273/NA); 10 mi E 
Rancho Melling (CIB 3367/MK872274/NA). 


Peromyscus truei montipinoris .—California; Kem Co., Rancheria Creek, east end Walker Basin (MVZ 197314/ 
MK871835/NA; MVZ197315/MK871836/NA; MVZ197316/MK871837/NA; MVZ197317/MK871838/NA); 
Temblor Range summit onHwy. 58 (MVZ 198606/MK871841/NA; MVZ198607/MK871842/NA; MVZ198608/ 
MK871843/NA; MVZ198610/MK871845/NA); 2 mi NNW Eagle Rest Peak, San Emigido Mts. (MVZ198616/ 
MK871846/NA); Los Angeles Co., Chatsworth Reservoir Park (TTU83290/AY376432/NA); 4.5 mi E (by 
road) Gorman (MVZ198394/MK871840/NA); 0.4 mi W Gorman (MVZ198392/MK871848/NA; MVZ198393/ 
MK871849/NA); San Luis Obispo Co., 13.3 mi NW (by road) New Cuyama (MVZ196792/MK871877/NA); 
Ventura Co., mouth Rose Valley (MVZ198613/MK871878/NA; MVZ198612/MK871880/NA). 

Peromyscus truei nevadensis .—Idaho; Cassia Co., City of Rocks National Reserve (UWBM79645/FJ800578/ 
NA; UWBM79646/FJ800579/NA; UWBM79674/FJ800582/NA); Nevada; White Pine Co., 5.3 km W Baker 
(BYU38274/MK871958; BYU38276/MK871960/NA); Utah; Beaver Co., San Francisco Mountains, 1.75 km 
N, 0.9 km E Frisco Peak summit (UMNH37585/MK871974/NA); Box Elder Co., 3.6 km N, 8.8 km W crystal 
Peak (BYU40195/MK871983/NA); Carbon Co., 3.3 km E, 2.5 kmN Castle Gate (BYU38484/MK871970/NA); 
Emery Co., 8.8 km S, 0.5 km W Window Blind Peak (BYU40714/MK871826/NA; BYU40706/MK871827/NA; 
BYU40704/MK872137/NA); 8.05 km N, 12.25 km W Huntington (BYU36684/MK872045/NA; BYU36685/ 
MK872046/NA; BYU36686/MK872047/NA; BYU38038/MK872077/NA); 7.95 kmN, 12.30 km W Hunting- 
ton (BYU36674/MK872065/NA); 12.5 km N, 22.7 km W Last Chance Benches (BYU38489/MK872104/NA); 
Garfield Co., Grosvenor Arch Day Use Area (BYU20180/MK872032/NA; BYU20181/MK872033/NA); Millard 
Co., Ferguson Desert Snake Pass Rd., 3.7 km E Shotgun Knoll (BYU24608/MK872107/NA). 


Peromyscus trueipreblei .—Oregon; Deschutes Co., 5 km W Tumalo (BYU21081/MK871963/NA). 

Peromyscus truei truei .—Arizona; Coconino Co., Navajo Nation (BYU41011/MK871823/NA; BYU41012/ 
MK871824/NA; BYU41013/MK871825); Apache Co., (Locality 13; TTU104427/AY376433/MN164793); Navajo 
Co., 3 mi S Woodruff (Locality 14; TTU78507/AF155412/MN164792); Colorado; Chaffee Co., 7.2 kmE, 3.65 km 
N Poncha Mountain (BYU37099/MK871883/NA); Mesa Co., W Grand Junction, Colorado National Monument 
(Locality 12; TTU55604/MN022896/MN164790); Fort Carson Co., Camp Red Devil (Locality 7; TTU100425/ 
MN022897/MN164791); Nevada; Clark Co., 1.20 kmN, 0.25 km W Willow Spring (BYU39546/MK871887/NA; 
BYU39550/MK871891/NA; BYU39552/MK871893/NA); 1.10 kmN, Willow Spring (BYU39542/MK871896/ 
NA); Willow Creek, 1.24 km N, 0.10 km W Willow Spring (UMNH40544/MK871899/NA; UMNH40545/ 
MK871900/NA); Spring Mountains, Telephone Canyon, 0.4 km S, 4.2 km E summit Fletcher Peak (UMNH40768/ 
MK871910/NA); Mineral Co., Wassuk Range, Cottonwood Creek, 2 mi S, 4 mi W Walker Lake (UMNH39130/ 
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MK871915/NA); Wassuk Range, Cottonwood Creek, 3 kmW summit of Mount Grant (UMNH39141/MK871918/ 
NA); Nye Co., Peavine Canyon, 0.92 km N, 0.27 km W mouth Horse Canyon (BYU33806/MK871922/NA); 
Pine Creek Canyon (BYU34412/MK871925/NA; BYU34413/MK871926/NA; BYU34415/MK871928/NA; 
BYU34417/MK871930/NA; BYU34418/MK871931/NA; BYU34434/MK871941/NA; BYU34442/MK871949/ 
NA); New Mexico; Catron Co., Quemado, Gila National Forest (Locality 11; TTU119068/MN022899/MN164796; 
TTU119069/MN022898/MN164795); Socorro Co., 32 mi S, 23.5 mi W Socorro (Locality 10; TTU36053/ 
AY376434/MN164794); Sandoval Co., 4 mi S, 3 mi WBemalillo (Locality 9; TTU36096/MN022900/MN164797); 
Harding Co., 9.0 mi SW Mills, Mills Canyon Campground, Kiowa National Grassland (Locality 6; TTU142741/ 
MN022901/MN164798; TTU142742/MN022902/MN164799; TTU142743/MN022903/MN164800); Oklahoma; 
Cimarron Co., Black Mesa Nature Preserve, (Locality 4; TTU136502/MN022904/MN164801; TTU120630/ 
MN022905/NA; TTU120631/MN022906/MN164802; TTU120632/MN022907/MN164803); Texas; Culberson 
Co., Upper Dog Ranger Station, Guadalupe Mountains National Park (Locality 8; TTU20584/MT670436/NA); 
Marcus Cabin, West Dog Canyon, 0.75 mi W, 6.38 mi N Guadalupe Peak, Guadalupe Mountains National Park 
(Locality X; TTU23544/MT670435/NA); Deaf Smith Co., 10 mi N, 35 mi W Hereford (Locality 5; TTU57016/ 
MT188565/NA); 11 mi S, 2 mi E Glenrio (Locality 5; TTU58045/MT188566/NA); Utah; Garfield Co., Escalante 
River Trailhead (BYU20114/MK872015/NA); 9.35 km E, 0.20 km N Mount Pennell (BYU35903/MK872115/ 
NA); 5.75 km E, 2.54 km S Steep Creek Bench (BYU35947/MK872186/NA); 11.00 km E, 0.90 km S Steep Creek 
Bench (BYU35977/MK872188/NA; B YU35978/MK872189/NA;B YU35986/MK872196/NA); 16.30 kmE, 5.85 
km S Steep Creek Bench (BYU36017/MK872199/NA); Wolverine Petrified Forest (BYU20316/MK872223/NA; 
BYU20323/MK872224/NA); 0.50 kmE, 1.30 kmN Wolverine Bench (BYU36047/MK872245/NA; BYU36055/ 
MK872252/NA); Grand Co., 0.25 km S, 0.30 km E Dewey Bridge, south of USA, Colorado River (UMNH38144/ 
MK872000/NA); 0.60 kmN, 0.90 km E Dewey Bridge, north of USA, Colorado River (UMNH38127/MK872001/ 
NA); Utah Bottoms, Dolores River (UMNH34767/MK872002/NA; UMNH34775/MK872003/NA; UMNH34776/ 
MK872004/NA; UMNH35706/MK872008/NA; UMNH35708/MK872010/NA); Grand Co., Rio Mesa field sta¬ 
tion, south side of Dolores River (UMNH38077/MK872006/NA); Kane Co., Buckskin Mountain (BYU23634/ 
MK871969/NA); Camp Flat (BYU23636/MK871977/NA); Dance Hall Rock (BYU20232/MK871989/NA); 
Devil’s Garden (BYU23591/MK871993/NA; BYU23574/MK871995/NA; BYU23577/MK871997/NA); 4.50 
km E, 5.30 km S Elephant Butte (BYU37451/MK872012/NA); 6.50 km E, 9.20 km S Flag Point (BYU37475/ 
MK872018/NA; BYU37476/MK872019/NA; BYU37479/MK872022/NA); Fourmile Bench (BYU23644; 
MK872027/NA); Kitchen Corral (BYU23656/MK872085); No Man’s Mesa (BYU23690/MK872122/NA); 
Smoky Hollow (BYU23731/MK872160/NA); San Juan Co., 0.8 km E, 3.6 kmN Navajo Mountain (UMNH43001/ 
MK871828/NA); 4.5 km E, 7.6 km S South Peak, Abajo Mountains (BYU38499/MK871968/NA); La Sal Mtns, 
tributary to Brumley Creek (UMNH31487/MK872096/NA; UMNH31406/MK872098/NA); Uintah Co., Rat 
Hole Canyon (BYU17166/MK872130/NA); Wayne Co., Aquarius Plateau, Carcass Creek, 3.0 km S, 0.5 km W 
Grover (UMNH31082/MK872040/NA); Torrey (UMNH36133/MK872210/NA). 
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